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PROBLEM WRITE UP 

1. Task Number: 59"05& Recoup II Computer 

2m Date Rec'd: 10 August 1959 

3* Submitted by: Mr. Mancini, Air fc Grnd Techniques Section, Suveyng & Good Br* 

4* Programmed by J L* A* Gambino 

5* Description of Problem: To determine the position of an unknown point with 

reference to the positions of two known points by 
making a continuous record of slant ranges to the 
three stations from an aircraft flying along an in- 
termediate path. Calculate the grid distances of 
this intermediate path using it as the base in the 
determination of the grid coordinates of ths unknown 
station* 

6. Mathematical Analysis: The trignometric formulae used in the position compu- 
tations were revised to include the concept from 
analytic geometry of the directed distance from a 
line to a point thereby making the formulae applicable 
to all cases* 

7* Numerical Analysis: Newton- Raphson iteration process was utilised to obtain 

approximate latitudes to be used in the scale factor 
confutations* The coordinates of the unknown can be eb_ 
talned to + 0.04 meters. 

8. Block Diagram: Attached 

9. Program Listing: Attached 

10. Operational Notes: 

a) Input and output via the electric typewriter* After reading in the tape, 
the location counter will be set at location 0001.0. Press the "Start 1" 
button on the right side of the console. The computer will now wait for 
input data* 

1. type A or B. Type A if the unknown station is above the known 
base line or type, F if it is below the base line. 

2* XXXXXXXX.XXX Ny (Nothing of the west point on ths bases* 

3. XXXXXX.XXX £g Basthg of the west point on the base. 

4. XXXXXXX.XXX (Ng) Northing of the last point on the base* 

5. XXXXXX.XXX (Eg.) Easting of the east point on the base. 

6. XXXX.XXX (fy) Elevation (in meters) of the west point on the base 

7. XXXX.XXX {#£) Elevation (in meters) of the root point on the base. 



6. XXXX.XXX 

9* XX.XXX 

10. XX.XXX 

11. XXXX.XXX 

12. XX.XXX 

13. XX.XXX 
14* XX.XXX 

15. XXXX.XXX 

16. XX.XXX 

17. XX.XX 
IS. XX. XX 

19. XX.XX 

20. XX.XXX 

21. XX.XXX 

22. XX.XX 

23. XXXX.XX 

24. XXXX.XX 

25. XX.XX 

26. XX.XX 

27. XXXX.XX 

28. XXX.XX 

29. XXXXX.XXX 

30. XXXXX.XXX 

31. XXXXX.XXX 

32. XXXXX.XXX 

33. XXXXX.XXX 

34. XXXXX.XXX 



hy) Elevation of the unknown station. 

TDBty°F) Dry bttlb temperature at vest end of baee. 

TWB^°F) Viet bmlb temperature at weet end of baee. 

ALfy ft.) Atiaete* reading (m feet) at west and of baee. 

Ety) Eccentricity at west end of baee* 

TDBg°P) Dryabulb temperature at eaat end of base. 

TWBg°F) Wet bulb temperature at eaat end of base. 

ALTg ft.) Altimeter reading in feet at eaat end of base. 

Efl) Eccentricity at east end of base. 

TOBy°F) Dry bulb at unknown point. 

WBy°F) Wet bulb at unknown point. 

ALTyft) Atlimeter reading in feet at unknown point. 

Ey) Eccentricity at unknown point. 

TDB°F) Dry bulb at west end of air baee. 

TWB°F) West bulb at weet end of air base. 

ALT) Altimeter, in feet at west end of air base. 

h) Altitude of aircraft in meters, at west end of air base. 

TDB°P) Dry bulb at east end of air base. 

TWBOF) Wet bulb at east end of air base. 

ALT) Altimeter, in feet, at east end of air base. 

h) Altitude of aircraft, in meters, at east end of air base. 

W-Air W) Distance, in meters, from West base paint to west end 
of air base. 

E- Air W) Distance from E ft st base point to west end of air base. 

V - Air W) Distance from unknown to west end of air 1 base. 

W - Air E) Distance from West base point to East end of air base. 

CS - Air E) Distance from East base point to east end of air base. 

(V - Air E) Distance from unknown point to east end of air base; 



Upon returning the carriage after enterintfHhis last distance , the 
computer goes into the compute mode* After approximately 5 jninutes, the 
position of the unknown will be typed out and the computer will wait for 
more distances to be entered* If there are no more distances, press the 
"Letters Shift" on the typewriter and the average of all the computed 
positions made on the unknown will be typed out. The computer will then 
set the location counter to 0001.0 and wait for a new problem. 

10b. 1. 1603 words or 3216 coiwaands are utilized by this program. 

2. Press the carriage return after typing in each piece of data. If an 
error is made while typing in a number, press the n J w key on the typewriter, 
press the start button on the console and then retype the number. If an error 
is detected after entering a nusfcer, it is advisable to read the tape back into 
the computer so that modified commands and stepping constants will be reset to 
their original values. 
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U. S. ARMY ENGINEER 
GEODESY, INTELLIGENCE AJID MAPPING RESEARCH AND DEVELOPMENT AOENCT 

suRVEnfrG and obodkst brahch 

1 lebruary 1961 
CALCULATION OF POSITION FROM DISTANCES MEASURED WITH 
AIRBORNE TELLUROMETER 



CONTINUOUS TRILATERATION METHOD 



1. General 



Electronic distance Measuring equipment, which has been given the 
distinguishing name ' Aerodlst ' , was developed in the Union of South 
Africa, for the U. S. Amy Itagineer Corpe (GIMRADA), for long range ext- 
ension of horizontal control. Essentially, the purpose of the equipment 
is to determine the position of an unknown point with reference to the 
positions of two known points. In the diagram below, A and B represent 
the two fixed ground stations with forward and back azimuths and C is the 
unknown point whose position is required. 




B 

A continuous record of slant ranges to the three stations A, B and C 
is made in an aircraft flying along an intermediate path, such as P Q. 
The aircraft is flown at as low an altiutde as possible, consistent with 
the maintenance of line of sight conditions. In flat terrain, the 



altitude would be of the order of 2000 foot for ranges of approximately 
50 miles. The height determinations of the aircraft and the unknown 
ground point are made barometrically. 

2. Computational procedure 

Of the several possible methods of calculation, that which sires accurate 
results most directly is the computation on the ellipsoid. 
The successive stages of computational procedure are: 

a. Selection from the continuous record of field data of at least six 
sets of simultaneously measured slant ranges. Six is considered a minimum, 
but provision for twelve should be made in the program. Six of these to 
air point f and six to point Q. For simple illustration in this problem 

for computational procedure, however, we shall consider only six} three would 
correspond to positions P.,, P~ and P- in the neighborhood of point P and 
three to positions Ql, Q , and Q« in the vicinity of point Q. Refer to 
Figure No. 1 for geometry of the problem. To insure optimum results, it is 
necessary that the figure and especially the triangle PQC be reasonably 
well conditioned. 

b. Computation of index of refraction using meterological 
data obtained at each point. 

c. Reduction of the slant ranges to ellipsoidal distances. 

d. Calculation of azimuths from known ground points A and B to air 
points P 1 , P 2 , P^ and a, Q 2 , Q^ projected to the eUpaoid and from Q^, Q^, Q~ 
to unknown ground point C. 

e. Computation of positions of air points and unknown ground point, 
which will be performed for each set of data with the final ground position 
taken as the mean of all C values. 



3. Reduction of slant ranges to spheroidal distances . 

The following procedure of reduction corrects the radio path *I/ to 
ellipsoid distance S for each measured line. Figure 2 contains the 
relationship of the variously computed distances to L. All lengths are 
expressed in meters, h, is the* elevation of ground point and hp the 

A\r PoioT 



elevation of aircraft (also in meters). 



Qrou/id foloi 




£ II ipso id , 



Figure 2. 



a » Reduction of L to P . 

The reduction of each radio path L to slant range D corrects for 
index of refraction (>-)) and is accomplished by application of the method 
contained in the Airborne Tellurometer Computations, Line Crossing Method 
(2 Feb 1961). The meteorological data will be measured and treated in the 
same manner as in the Line Crossing Problem. 

b. Reduction of D to 15 . 

Slant ranges are next reduced to slant chord distances D for each 
line AO, AP 1 , BO, BP-, Q.O., P.C., etc., by the formula 

D-D - D 3 /&r 2 (1) 

where r - 26.0 x 10 meters (radius of radio wave). 
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c. Reduction of D to K . 

1. In performing the reduction of D to ellipsoidal chord distances 
K, approximate asimuths of all measured lines and geographic coordinates of 
points P., P 2 , Q., Q-, etc., must first be determined. By giving the various 
distances their respective asimuths prior to further reduction makes possible 
a refinement of the distance since in many cases these measured lengths are 
likely to be in excess of 100 miles. Formulas contained on page 1 of Problem 
Mo. 8, Trilateration Adjustment (3 April 1956) of the Engineers Problems can 
be used for computing internal angles of the triangles of the quadrilateral. 
The asimuths of lines AQ., AQ~, BF., BP^, etc., need only be computed. 
Asimuths may be added according to the convention set forth in the Line 
Crossing Method. 

2. Hext, the approximate geographic coordinates of the air points P., , P-, 
Ql, Q_' •*©.* are required. These coordinates can be determined by employing 
Helmert's Position computation which is a part of this problem and utilises 

as input data initial coordinates of A(orB) and asimuth to the airpoint Q(orP). 

3. following these preliminary computations! the reduction of D to K is 
relatively simple. By application of the following formula utilising asimuth 
and position data thus determined, the ellipsoidal chord distance K is obtained. 



i 



2 -fa,-**)* ' ( 2) 



i P CT 

where D, h. , and h* ane known, 

and 

p - pv/(v cos a + p sin c) - radius of curvature. 

also p ■ a (l-e ;/fl-e sin ) 

mm 

a s l /z 
v ■ a/(l-e sin M ) 

further a - 6378206 m - semi-major axis of ellipsoid. 

t - latitude of starting point of line, 
m 

a - asimuth from A to a, B to P., etc., from north 
(previously computed). 



d. Reduction of I to 8, 

Finally, to obtain ellipsoidal arc distance S, we use the relation 

8 " * + il * il£ (3) 

*K* °uoC 

Where K and p are defined above, 

e. Adjustment of triangles. 

Each triangle ABP-., ABP^, ABO, ABQ*, etc., ia now ready for 
final adjustment. By utilising the formulae 

cosB, - cos H - cos SP cos SI (U) 

3 Sin BP Sin Si! 



COSA • cos SP - COB 51 COB PI (5) 

Sin SI Sin PI 



cosP.,- COS 15 - COS IP COB SP (6) 

Sin IP Sin SP 



; > 

compute the angles of the two triangles the sum of each should - IdO + ( , 

where f- spherical excess to be computed according to eq. (lU) in Line Crossing 
Method, with t - mean of vertices of each triangle. 

Since this condition rarely occurs, there will normally be an error of closure E 
of the triangle. The error of closure is to be eliminated according to the pro- 
cedure in Line Crossing Method, page 7. 

f. Final asimuth . 

Final asimuth of lines BP-, BP 2 , AQ^, AQ 2 , etc. can be determined by add- 
ing (or subtracting) the angles A- , A, , B, , B- , etc., to the asimuth of 
line AB (known) according to the convention on page 7 of the Line Crossing 
Method. 
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g. Computation of line PQ . 

1. The final geographic coordinates of points ^0%$ P 2 » ®2> ** c »> 
are to be conputed from the final azimuths above and the ellipsoidal lengths 
'S* previously computed. Helmert's Position Computation is to be used. 

To compute the length P«iQj> *V^2* e ^°*> * 0O P7 °^ So^aa©** inverse 
formula (with accompanying example) is also furnished and will be used for the 
computation. This formula supplies, equations for computation of forward and back 
azimuth 4 also. 

h. Computation of unknown point C . 

1. From the computed lengths and atimuths of line PQ, the approximate 
azimuths of lines QkC-, QC-, J^CL* ?o®2' etc *> * ra computed along with coordinates 
for C*, C 2 , etc. Internal angles of the triangle PjQiCj, ^W %te *> art 
adjusted in the same manner as above for the triangles of the quadrilateral, 
except this time lines P-tQk* p o^2' *tc.j are considered as the base line for 
each set of data instead of AB as before. 

2. Final azimuths &C,, Q 2 C 2' • tc «* ar * determined as above in f. 

3. The geographic coordinates of point C., C-, etc., from Q., Q^, etc., only 
are computed and finally, the latitudes and longitudes of the various values of 
C are meaned to determine the most probable value. 

i. Synopsis . 

Since the problem is set fourth in some detail, a quick review of 
computational procedure for single value of C is felt desirable. This review 
is supplied as an aid to the programmer in hopes of simplifying the problem. 

From the given data of stations A and B, coordinates of air points P, and 
Q. (reduced to the ellipsoid) are oomputed over lines BP, andiQ. as a function 
of azimuth and distance. An inverse of line P, <1 is performed, and using this 
line as a base with azimuth Q. C. and the related distance, geographic coordinates 
for C. are computed. This procedure is carried out for each set of measurements 

resulting in a unique value for C. The C's are then meaned as described above 
to yield the most probable value for the unknown point. 

j. Program flexibility . 

The program will be developed for optional input and required output 

of the following basic geodetic data. This problem operates on the ellipsoid 
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( geodetic); and should input data of tha starting points ba on tha UTH 
plain or In mils, than It will be required that these data be converted 
to geodetic quantities ria subroutines (integral part of the program). 

1. Input . 

a. Geodetic 

(1). Geographio coordinates ( ° ' * ) 
(2). Azimuths ( ° ' * ) 

b. UTM . 

(1). Coordinates (meters) and zone, converted by subroutine to 
geographic coordinates. 

(2). Azimuth ( ° * * on mils), converted by subroutine to 
geodetic azimuths. 

2. Output . 

a. Geodetic 

(1). Geographic coordinates ( ° * * ) 

(2). Azimuth ( ° > * ) 

b. UTM . 

(1). Coordinates (meters) and zone. 
(2). Azimuth ( ° ' * and mils ). 

k. Utilise formulae of Engineer problem no. 12 (UTM to geographica) to 
convert UTM coordinate input to geographies. From other existing formulae convert 
UTM sexagesimal or mil azimuth to geodetic. 

1. In output utilize Engineer survey problem 11 formulae for Geographic 
•to UTM coordinates. From convergence formulae, convert geodetic azimuth to 
UTM sexagesimal and mil azimuths. 

m. Entire program will contain as integral part all subroutines. 

n. Program will ba designed to work only north of the equator and west 
of Greenwich. 

o. Program shall ba developed so that following entry of Baseline data, new 
Field data can be entered for computing other C's. 

Motet Refer to Format as guidance. 



